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We propose a simple approach to study the conductance through an array of A'^ interacting 
quantum dots, weakly coupled to metallic leads. Using a mapping to an efTective site which describes 
the low-lying excitations and a slave-boson representation in the saddle-point approximation, we 
calculated the conductance through the system. Explicit results are presented for A'^ = 1 and 
A'^ = 3: a linear array and an isosceles triangle. For A'^ = 1 in the Kondo limit, the results are in 
very good agreement with previous results obtained with numerical renormalization group (NRG). 
In the case of the linear trimer for odd A'', when the parameters are such that electron- hole symmetry 
is induced, we obtain perfect conductance Go ~ 2e^ /h. The validity of the approach is discussed in 
detail. 
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I. INTRODUCTION 

Potential technological applications in electronic de- 
vices has led to intense research in electronic transport 
in nanoscale systems. In particular, a quantum dot (QD) 
acts as a single electron transistor and in addition, the 
many-body physics of the "single impurity" Kondo effect 
has been clearly observed in transport experiments. 
More recently a system of two QD's has been studied, 
which provides an experimental realization of two Kondo 
"impurities" interacting with a metallic host.'' Also, lin- 
ear arrays of 15 QD's have been fabricated to investigate 
its electronic properties,^ and the conductance through 
a linear trimer of QD's has been investigated in the 
Coulomb blockade regime.^ Systems of coupled QD's are 
also of interest because of their possible application in 
quantum computation.'' 

Theoretically, the conductance through one QD at 
equilibrium has been clarified by different studies of 
the impurity Anderson model, in particular those using 
the highly accurate Wilson's numerical renormalization 
group (NRG).®"^" For moderate values of the Coulomb 
repulsion U, perturbation theory in U (PTU) provided 
good results. ^^^^^ In particular, an approach that in- 
terpolates between the second-order results for the self 
energy and the atomic limit, '^^"'^■^ provided accurate re- 
sults, as has been shown by direct comparison with ex- 
act diagonalization in small rings. Renormalized'^^ and 
intcrpolativc^^ PTU have also been used in the non- 
equilibrium case. At zero temperature, the method of ex- 
act diagonalization plus embedding (EDE) was success- 
fully used.^*"^° This method is based on solving exactly 
some part of the system, which includes all interactions, 
and embedding it in the rest of the system. A detailed 
study of the approximation and its range of validity was 
given by Chiappe and Verges. A brief description of the 
PTU, EDE and slave boson approximation is contained 
in Ref. 22, where the approaches are applied to describe 
scanning tunneling spectroscopy of different systems with 
magnetic impurities on Cu or noble metal (111) surfaces. 

Transport through a system of two QD's, one of them 



coupled to two conducting leads has been studied using 
NRG.^^'^'' More recently NRG rcsuhs for three^^ and 
four^^ QD's on a line have been reported. These two 
works deduce the phase shifts for even and odd pari- 
ties from the effective non-interacting Hamiltonian that 
describes the strong coupling fixed point. For Fermi liq- 
uids with inversion symmetry, the conductance can be 
expressed in terms of these two phase shifts. In spite 
of these studies, for several QD's or systems with low 
symmetry, application of the NRG becomes impractical 
and one has to resort to other approximate techniques. 
The EDE has been applied recently to several problems 
involving more than one QD.^^"'^". In particular, Biisser 
et al. studied the conductance through a linear array of 
QD's.^^ Remarkably, they find that for an odd number N 
of QD's except N = 1, the conductance G vanishes in the 
electron-hole symmetric case (EHSC). For = 1, G as a 
function of gate voltage Vg reaches the maximum at the 
EHSC with the ideal value Go = 2e'^/h in agreement with 
previous studies. ^"'^^'^^■-'^'^^ However, for = 3, 5, ... the 
peak is split in two by a deep minimum. This result is 
in contradiction with previous results using PTU which 
predict ideal conductance G = Go at the EHSC.'^^ There- 
fore, Biisser et al. conclude that other many-body tech- 
niques should be used to elucidate the issue. 

The purpose of the present work is to present an ana- 
lytical approach to calculate the conductance through an 
interacting region, weakly coupled to conducting leads. 
As an application, we study the above mentioned con- 
troversy. In particular, we consider an interacting region 
composed of a linear array of an odd number N QD's 
(extension to even N and other geometry is straight- 
forward), connected to two non- interacting leads by the 
same hopping V. We also study the conductance for the 
case in which the interacting region consists of an isosce- 
les triangle of QD's. To simplify the discussion when the 
energy of states with different number of particles are 
compared, we set the origin of one-particle energies at 
the Fermi level ep = 0. In general, if the ground state of 
the interacting part is quasi degenerate between a non- 
magnetic singlet and a magnetic doublet, and V is small 



compared with the difference in energy with other states, 
neglecting the latter the problem can be mapped into 
a one-impurity Anderson model. As stated above, this 
model is well understood. This mapping has been used 
to calculate the conductance through a ring described by 
the ionic Hubbard model, which should detect a topolog- 
ical phase transition. '^^ In the Kondo regime, where the 
magnetic ground state lies well below the other states, it 
might be necessary to retain other states with even total 
spin for an accurate description.'^^. As in these problems, 
we assume that the ground state of the interacting part 
is a doublet for an odd number of particles. In particu- 
lar, the Hubbard model used to represent a linear array 
of N QD's has a doublet ground state in the EHSC for 
odd N. We retain this doublet and all singlet states 
which hybridize with it to map the problem into a gen- 
eralized one-impurity Anderson model. The validity of 
the mapping is discussed in detail in section IV. The re- 
sulting model is solved using a slave-boson representation 
in the saddle-point approximation. Explicit calculations 
are presented for A'^ = 1, a linear trimer and an isosceles 
triangle. The main result is independent of odd N for 
a linear array of QD's: we obtain ideal conductance in 
the EHSC. This agrees with recent results obtained us- 
ing alternative techniques. ^^'^^'"^^ The isosceles trimer is 
of interest because of the effects of magnetic frustration 
and peculiar electronic structure. In this case the 
electron-hole symmetry is lost and this fact is reflected 
in an asymmetric lineshapc for the conductance. Despite 
this fact, the unitary limit is reached when Vg favours a 
magnetic doublet ground-state. 

The paper is organized as follows. In Sec. II, we 
describe the model, the mapping to one effective site, 
and the slave-boson formalism for the resulting effective 
model. The results for the conductance are contained in 
Sec. HI. Sec. IV summarizes our results and discusses 
our findings in relation to previous works. Some details 
of the calculations are left to the appendices. 



II. MODEL AND APPROXIMATIONS 



Lead Interacting Dots Lead 
t_t„V^t>-^t>-NV„ t_ t 




Figure 1: Scheme of the studied system. 
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ffd describes the central region containing the dots, each 
with an on-site repulsion [/, 
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Finally the term 77mix, that couples the interacting part 
of the system with the conducting leads has the form 



H.C. 



(4) 



We take the Fermi energy ep=Q, in such a way that the 
leads are half filled. If in addition Vg is such that Ed = 
—U/2, the Hamiltonian (for a linear array of QD's) is 
invariant under an electron-hole transformation ^ — > 
(— l)'ci.o-. In addition, for any value of Vg, the system is 
invariant under inversion cj^ cl^ ^ (for odd N). 



A. Model 

For the sake of clarity in the presentation, we consider 
a linear array of an odd number N of QD's coupled to 
the left and to the right to metallic leads, and take = 3 
unless otherwise stated. The changes for any other array 
of QD's are obvious. The system is represented in Fig. 
1. The dots are equivalent and their on-site energy Ed = 
—eVg is controlled by the gate potential Vg. The leads 
are described by non-interacting half-infinite chains. 

The Hamiltonian is 

H = i/i + + i/mix. (1) 



B. Mapping of the interacting region to an 
effective site 

We assume that the state of lowest energy of the in- 
teracting part of the system for a certain odd number of 
particles of interest is a doublet. This is certainly true 
for the ground state of Eq. (3) in the EHSC. i/mix mixes 
this doublet with other states with A — 1 and A -I- 1 
particles. For small enough V, -ffmix can be eliminated 
trough a canonical transformation, and the model can 
be mapped to a single impurity Kondo system, '^^ which 
is well understood. However, as Vg varies, the state of 
minimum energy in the sector of either A — 1 or A -I- 1 
particles becomes quasi degenerate with the ground state 



2 



and finally crosses it, rendering the canonical transforma- 
tion invalid. Near the quasi degeneracy the model can be 
mapped into a one-impurity Anderson model. '^^ Here we 
generalize previous approaches'^^''^^ and map the system 
into a generalized Anderson model, retaining the lowest 
doublet in the subspace of TV particles, and all singlet 
states in the subspaces with iV ± 1 (see Fig. 2). For 
simplicity we neglect the triplet states, but as we will 
show this does not modify our conclusions. The neglect 
of excited doublets is the most serious and its effect is 
discussed in Sec. IV. 

The eigenstates of are classified according to its 
parity (z/ = ±), total spin S and its projection 5"^. We 

will denote these eigenstates as: IV'j^^): where n is the 
number of particles, j orders the states according to in- 
creasing energy and p is a short label for the other quan- 
tum numbers {u, S and Sz)- The corresponding energies 



are denoted by E[ 



n=2 subspace 
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with the following constraints 
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n=3 subspace 



n=4 subspace 



We point out here that for a linear trimer, the ground 
state \^pQ _ g.) of the interacting part is odd under in- 
version. This implies that one has to choose one of the 
operators f^, even and the other odd under inversion. 
Otherwise the invariance under inversion would be vio- 
lated. 

In this representation, the Hamiltonian can be written 

as 




Figure 2: Scheme of the low-lying states of for N = 3. The 
dashed lines denote states truncated in the approximation. Dashed 
arrows represent processes of order V*, which are neglected. Only 
processes of order V'^ are retained. 

To represent the truncated Hamiltonian, we will use 
a slave boson formalism. This is convenient in our case 
in which, in addition to the truncation, we perform a 
saddle-point approximation to the slave boson formalism. 
However, the slave boson representation is independent 
of the method of solution chosen, and has been used for 
example in exact diagonalization of finite clusters includ- 
ing triplet states. '^^ 



C. Formalism of slave bosons 

As stated before, for the sake of clarity we will refer 
to the case of the linear chain. We represent the many- 
body states of i?d using a slave boson representation that 
respects electron-hole symmetry. This symmetry is im- 
portant for the description of the conductance through a 
linear array, as discussed in the introduction. Basically, 
this representation can be described as a generalization 
of that used by Kotliar and Ruckenstein for the Hubbard 
model^° to include several "empty" e and "doubly oc- 
cupied" d bosons. Specifically, we introduce the boson 
operators ej^^,dj^i, and and the fermionic ones f^, to 
map the states of the effective site describing the inter- 
acting part as 

'.(2)\ 



1^ 



eL|0), 



H — Hi + H[ + H„ 



He, 



(8) 



where Hy is the same as before (Eq. (2)). Hf describes 
the energy of the fermions at the effective site 



(9) 



where S^^) ^ £^(3^ ^ jg ^-^^ 

energy of the lowest doublet 
in the subspace of three particles (we keep = 3 in this 
Section for simplicity). Similarly for the bosons at the 
effective site 



2^(3) 



(10) 



The energy of the boson Sg. has been chosen at zero. Note 
that the above choice ensures the correct energy for each 
state of the effective site: 

PoHPo eljO) = Efl e],jO), 

Pi if Pi /t4|o) = /t4|o), 



{J^^-2P^) + 2P^ 
x4„/|/||0). 



(11) 



where we have used the projectors P,; over the subspace 
with i fermions, orn = i-|-2 (n = i — 1 + A^in general) 
particles in the interacting region. 

Defining for convenience electron operators with well 
defined parity under inversion 



V2 



(12) 



the term -ffmix in this representation becomes 
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V2V { J2 
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(13) 
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where 



(14) 

= (d4,(...,),jC^>' (15) 

where is the parity of the ground-state for N particles. 
For the Hnear array, ofj-.i/ and (3j^i^takes the same values 

I 



',(2)\ 



as a consequence of the electron- hole symmetry for a par- 
ticular Vg (as mentioned above) and the fact that both 
matrix elements turn out to be independent of Vg. 

The constraints are incorporated in the Hamiltonian 
introducing Lagrange multipliers A' and Ao-, correspond- 
ing to the Eqs. (6) and Eq. (7) 



(16) 



In the functional integral formalism, the fcrmions can be integrated out as in Hubbard model, and the partition 
function becomes 



Z = j l[[Ve,,,]Y[[Vd,,,]l[[Vs„]dX,dX' cxp - j^r S{i 

where the effective action is 

Sir) = E + E^^l + A') e,.J\ + Y ['I + E(^^ + A' - A,) s, 



E 



dl^ {dr + Efl - 2ii;(3) + A' - ^ A,) d,, , 



- A' + ^Tr InM,^,, 



and the non-zero matrix elements My,„ (for iV = 3) are 



= { (a, - M + - A,) 



if i = ? 0, 

if (i = 0;/ = 2) or (i = 2;/ = 0) 



r 



(17) 



(18) 



where 



K,. - V2F 



(19) 



Up to now, the only approximation made was the trun- 
cation for iV > 1 of the Hilbert space, to a set of rele- 
vant low-energy states of the interacting part H^- For 
= 1, the formalism introduced in this subsection is 
just a change of representation of the original Hamilto- 
nian. 



To solve the problem we used the saddle-point approx- 
imation to evaluate the partition function. The bosonic 
fields are replaced by real constant numbers that mini- 
mize the action. The details are left for appendix A. 



III. RESULTS FOR THE CONDUCTANCE 

In the slave-boson mean-field approximation that we 
are using, the many-body problem is reduced to a non- 
interacting one with one "impurity" and only one chan- 
nel of conduction, with renormalized hybridization (see 
Eq. (19)). Then we can use the two-terminal Landauer 
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equation'*^ obtaining 



dco I \tiu;)\\ (20) 



where /(w) is the Fermi function and t{Lu) is the trans- 
mittancc through the effective site 



(21) 



Gffa{i^) is the retarded Green's function of the fermion / 
with spin a (see appendix A), and V±2 = {V+ ± V-) /\/2 
is the couphng between the effective site and (for = 3) 
the site with subscript ±2. At zero temperature from the 
above equations one gets 



G 



2e2 



[2npo{V^~V^)]'\Gff.\l^o- (22) 



A. Conductance for N = 1 

We first discuss the case of the conductance through 
one QD, and we compare the results with those obtained 
with NRG^ for different vahies of A, which is the half 
width at half maximum of the resonant level in the non- 
interacting case. In our problem A = 27TpoV^, where 
po = l/(27ri) is the unperturbed density of states at the 
Fermi energy for one spin and one of the leads, and the 
factor 2 adds the contribution of both leads. 

The results for one dot are shown in Fig. 3. The ideal 
conductance for Ed = —U/2 indicates the formation of a 
virtual bound state (Kondo effect) at T = 0. As can be 
seen from the inversion symmetry of the figure around 
Ed = —U/2, the solutions of the saddle-point equations 
reflect the electron-hole symmetry of the formalism. A 
very good agreement between our results and those of 
NRG for A/ttU = 0.01 (strongly-interacting case) is ob- 
served. This is in part due to the fact that in the EHSC, 
where analytical results can be obtained (see appendix 
B), the formalism reproduces the correct exponential de- 
pendence of the width of the resonance in the spectral 
density of states S, which roughly coincides with the 
Kondo temperature Tk- However, the shape of the res- 
onance is Lorentzian in the mean field approximation, 
while the correct result in the Kondo regime is^^ 



Pda{^) = —T-Re^/Juj+lSjJiS, 
ttA 



(23) 



where S = 1.55T/f .i" 

The good agreement between the slave-boson mean- 
field approximation and NRG deteriorates as smaller val- 
ues of U are considered and is lost in the non-interacting 
limit. Kotliar and Ruckenstein proposed a remedy to 
this problem, replacing some operators entering i/mix in 
Eq.(13) by another ones (containing some suitable cho- 
sen roots) which coincide with the previous ones when the 



constraints are imposed exactly on each site, and at the 
same time reproduce in the saddle-point approximation, 
the correct results for the noninteracting case.*° Unfor- 
tunately, we find that this method does not work in the 
whole range of A/nU, giving an overestimated Kondo 
temperature in the strongly- interacting limit. In con- 
trast, as mentioned above, the formalism so far presented 
gives the correct dependence of Tr- on the parameters in 
this limit (see appendix B). 
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Figure 3: Conductance vs E^/U ratio, for different values of 
A/ttU . Lines indicate the results of the present slave boson (SB) 
approach, and the square points have been obtained with NRG.* 

Since we are mainly interested in the limit of strong 
correlations, the results shown in Fig. 3 are encour- 
aging. For small or moderate values of U one can use 
PTU.^^'^'^'^'*'^^ In spite of this good agreement, the con- 
ductance falls to zero abruptly at some value of Ed- This 
is an artifact of the saddle-point approximation, reported 
before in other problems, for example in the con- 
ductance through rings of quantum dots.'^^ Beyond some 
value of the gate voltage Vg , the impurity decouples from 
the rest of the chain. 

The decrease in the conductance out of the EHSC is 
due to the reduction of the density of states at the Fermi 
level, as the system approaches the intermediate valence 
regime. Technically, in the saddle-point approximation 
as Ed increases, the magnitude of boson e also increases 
(the charge in the dot decreases). This causes the mag- 
nitude of the boson s to decrease and the resonance in 
/-electron density of states must be shifted to higher fre- 
quencies in order to decrease the n j occupation, accord- 
ing to the constraint Eq. (7). This causes a reduction in 
the value of IG/zo-l at the Fermi level, and a reduction 
in the conductance (see Eq. (20)). In a similar fash- 
ion, when the dot is charged, the value of the boson d 
increases, the value of n/ approaches to 1, and the res- 
onance in the density of states must be shifted to lower 
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frequencies, producing an analogous decrease in the con- 
ductance. 

For the simple case of conduction through one dot, the 
expression of the conductance can be written as^^'**^ 

G = ^Jda; (^-^^ ttA p,,{u;) , (24) 

At very low temperatures (T <^ Tk), the system is a 
Fermi liquid. Therefore, the Friedel's sum rule (FSR)^'* 
must hold. For the case in which the hybridization V 
and the unperturbed density of conduction states po do 
not depend on energy, the FSR takes the form 

sin^(7r(nda)) , , 

Pda(O) = T , 25) 

ttA 

giving a conductance equal to 

G = ^sin2(7r(nd,)), (26) 

where (n^cr) is the total charge per spin in the dot. 

It can be shown (see appendix B) that the saddle-point 
approximation verifies the FSR, since the real problem 
is mapped into an effective non-interacting one whose 
ground state is a Fermi sea. 

B. Conductance for A'^ = 3. Linear trimer 

The results for conductance through three QD's using 
the saddle-point approximation described in appendix A, 
are shown in Fig. 4. Results for one QD for the same 
parameters are also shown for comparison. 




Figure 4: Conductance vs Vg/U ratio, for different values of A/ttU. 
Innermost curves correspond to the conductance through 3 QD's. 
The 1 QD conductance is also shown for comparison. Parameters 
are U = t, t' = 0.5 t and the hybridization V is 0.25 t, 0.3 t, 0.4 t 
and 0.5 t for increasing values of A/nU. 

On a qualitative level, the results for three QD's are 
similar to those of one QD near the EHSC. In particular, 
the conductance equals the ideal one Go in the EHSC, 



decreases as the amount of charge fluctuations increases, 
and finally becomes very small when the ground state 
of the interacting region becomes non-magnetic. The re- 
duction of the region in which the conductance is near 
to the ideal one is reduced in the three-dot case, and one 
expects further reduction as the odd number of dots N 
increases. This is due to the fact that as N increases, 
the energy of singlets with ± 1 particles become nearer 
to that of the ground-state doublet with A^ particles. As 
a consequence, the interval of gate voltage Vg between 
both intermediate valence regimes (when the energy of 
one singlet coincides with that of the doublet) is reduced. 

When Vg is moved further from the EHSC, an impor- 
tant difference between A^ = 1 and odd A^ > 1 is that in 
the latter case, magnetic states with A^ ± 2 particles be- 
come accessible, a new Kondo effect takes place, and an- 
other plateau with near ideal conductivity appears. For 
N > ^ this pattern with regions of nearly ideal or zero 
conductivity can be repeated several times. In general, 
including even A^, there should be plateaus when the gate 
voltage favors an odd number of particles. In particular, 
for N = 2 the conductance in a similar system has been 
calculated using NRG and displays two plateaus, corre- 
sponding to 1 and 3 particles in the interacting region. 
While we have limited our calculations to odd A^ and 
gate voltages such that either A^ or A" ± 1 particles are 
favored, it is straightforward to extend the formalism to 
other cases, as long as the hybridization V of the inter- 
acting region with the rest of the system is small enough. 
The allowed magnitude of V for the validity of the for- 
malism is discussed in the next section. 

We have verified numerically that Eq. (26) also holds 
for three dots. The formalism is identical for other odd 
N. Therefore, near the EHSC the situation for different 
odd A^ are qualitatively similar in that they show a local 
Fermi liquid behavior, ideal conductance at the EHSC, 
and the expected electron-hole symmetry. This is at vari- 
ance with the results of Biisser et al.^^ who find vanishing 
conductance at the EHSC, but are consistent with those 
of Oguri,^'* and recent studies for A^ = 3,^^''^^ who also 
obtain ideal conductance at the EHSC. Note that the pa- 
rameters with V ~ 0.3 t used in Fig. 4 correspond to one 
set of parameters used in Fig. 2 (f ) of Ref. 28. 



C. Conductance for A'^ = 3. Isosceles triangle 

In this subsection, we consider a system with the ge- 
ommetry schematically depicted if Fig. 5. Two regimes 
in parameter space were analized, namely t" < t' and 
t" > t', where t" is the new hopping term connecting dot 
1 and -1. These regimes correspond to the conductance 
through the two possible "isosceles triangles" . In the case 
t" = t' (equilateral triangle), the subspace with n = 3 has 
a doubly degenerate ground-state, with states belonging 
to the even and odd subspace respectively. This can be 
understood by recalling that the symmetry group of the 
equilateral triangle in two dimensions {C^y ) has one two- 
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dimensional irreducible representation. As it will be dis- 
cussed in detail in Section IV, the present formalism is 
not valid at this point of parameter space, since we as- 
sumed a non-degenerate ground-state. Since experimen- 
tally is very hard to reach such a regime, we believe that 
this is not a serious constraint of the method. Starting 
with the linear chain {t" = 0), we recall that the ground- 
state of the n = 3 subspace is an odd S — 1/2 doublet. 
Switching on adiabatically the hopping t" , the energy of 
the excited even doublet begins to fall, and after crossing 
the degeneracy point t" = t', the even doublet takes over 
and the formalism is valid again. 

Besides these features, we expect the development of a 
plateau in the conductance at the ideal value in a certain 
region of Ed/U, due to the fact that the problem still can 
be mapped into an effective 5 = 1/2 Anderson impurity 
coupled to conducting leads. 
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Figure 5: Scheme of the triangular arrangement. 



Figure 6: Results for the conductance through the isosceles tri- 
angle for the case t" < t' . Parameters are t' = 0.5 t, U = t and 
V = 0.45 t. 



evolution of the points E^^ and E^^ as a function of t" 
(see Fig. 7), the evolution of the conductance curves in 
Fig. 6 can be explained easily: the conductance is related 
to the low-energy physics of the isolated trimer. For ex- 
ample, for t" = {t" = 0.4 t) , the first value of Ed/U for 
which G = Go/2 in Fig. 6 is near —0.8 (—1.2), the value 
for which the number of particles in the ground-state of 
the isolated trimer changes from 4 to 3 (see Fig. 7). As 
it is already clear from Fig. 4, the width of the crossover 
regions increases with V, since the effective hybridiza- 
tion Vj/^o- is proportional to it and therefore the widths of 
intermediate valence regimes increase accordingly. 



In Fig. 6 we show the results for the conductance 
through the isosceles triangle. For simplicity we have 
set parameters V = 0.45 t, t' = 0.5 t and U = t, so that 
the first result {t" = 0) is one of the curves of Fig. 4. 
As t" is increased, the curves are approximately rigidly 
shifted to lower values of Ed/U. This feature can be cor- 
related with properties of the isolated trimer. In order to 
gain physical insight, we define the following functions of 
the isolated trimer: 



A23iEd,t",t\U) = El'\Ed,t",t',U) ~ El^\Ed,t",t',U), 
A^siEd,t",t\U) = El^\Ed,t",t\U)-El^\Ed,t",t',U), 



which are the differences in energy between ground- 
states belonging to subspaces with different particle num- 
ber. Fixing t",t' and U, we define E^''-^ and E"^'^ as 
the values of Ed which render A23{Ed, t" , t' , U) ~ and 
A43{Ed,t" ,t' ,U) = respectively. These values corre- 
spond to an intermediate valence regime of the efective 
Anderson model, and are related to the crossover regions 
in the dependence of the conductance G with gate volt- 
age, where G falls down from its ideal value at the plateau 
to zero. This is due to charge fluctuations and the conse- 
quent dissapearance of the Kondo effect. Following the 




Figure 7: Evolution of the degeneracy points as a function of t" 
for the isolated trimer. Parameters t and U are those of Fig. 6. 

The results for t" > t' are shown in Fig. 8. The asym- 
metry in the curves is more evident, and again the plot 
can be understood by following the evolution of E^'^ and 
E^^ (see Fig. 9). When t" > t', the low-energy physics 
of the system resembles that of a singlet formed by the 
two QD's connected by t" plus an additional QD on the 



7 
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Figure 8: Results for the conductance through the isosceles tri- 
angle for the case t" > t' . Parameters are the same than Fig. 
6. 



-0.4 



LU 




P -1 .0 



Figure 9: Same than Fig. 7 for the case t" > t' . The different 
evolution of the degeneracy points as a function of t" explains the 
asymmetry of the conductance in Fig. 8. 



top. In this situation, there is a considerable decrease 
in the effective hybridization V,^.a- (see Eq. 19). This is 
directly reflected in that the Kondo temperature is much 
lower than before, rendering the convergence of the nu- 
merical method difficult. For that reason, in those cal- 
culations we have set V ~ 1.0 i in order to have a larger 
Kondo temperature. The validity of the formalism is only 
shghtly affected, despite the fact that V is comparable to 
other energy scales (t, U). The main reason for this is the 
high values of t" , which ensure a large difference between 
the ground state energy and that of the excited odd dou- 
blet for n — 3. Another reason is the small effective 
hybridization. Other details on the validity of the ap- 
proximation will be given in the next section. As a final 
comment, we mention that the results obtained for these 
geommetries are consistent with those obtained with the 
FSR (Eq. (26)) for the effective Anderson model. 
IV. DISCUSSION 



agonalizing exactly the part of the Hamiltonian describ- 
ing the interacting region and retaining the most relevant 
states which describe the low-energy charge excitations. 
As discussed below, the approach is valid when the hy- 
bridization V of the interacting region with the rest of the 
system is small. The advantage of this mapping is that 
the physics of one impurity is rather well understood, 
and several good approximations can be used. Recently 
a mapping to an effective Kondo model with one "impu- 
rity" has been successfully used to interpret the depen- 
dence of geometry of spectroscopic experiments for a Cr 
trimcr on Au(lll),'^^ and to calculate the conductance 
through a dimer connected at one site to two conducting 
leads. In the latter case, NRG results are available and 
are in excellent agreement those of the effective model 
for small or moderate y.^^ 

For small V, the effective site approach can also be 
used when the the interacting region is side coupled to 
a conducting lead and the current does not necessarily 
flow through the interacting region. An example of this 
system is a trimer coupled to one site of a conducting 
lead.^^ In this case, the model is mapped to an effective 
site side coupled to a conducting lead,^^'^^'^^'^^ instead 
of embedded in the lead. 

We have introduced a slave boson representation for 
the description of the effective Hamiltonian, and used the 
saddle-point approximation to calculate the conductance 
through the effective site. The results obtained in the 
strongly-interacting regime for N = 1 were found to be in 
very good agreement with NRG calculations. For a linear 
array of an odd number of QD's, our results near the 
EHSC are qualitatively similar. In particular, we obtain 
perfect conductance at the EHSC. Since this result might 
seem controversial at present, we discuss in detail the 
approximations involved in our treatment, particularly 
in relation to this result. 

The crucial approximation is to retain only the low- 
est doublet of the magnetic configuration, neglecting the 
excited ones (see Fig. 2). This approximation scheme 
is only valid as far as the excited levels have a negligi- 
ble contribution to the many-body ground state. The 
hybridization V induces a second-order matrix element 
Af°5>e between the ground state of the magnetic config- 
uration in the interacting region and excited states for 
the same number of particles. In particular, for the lin- 
ear trimcr, while the ground state for three particles is 
an odd doublet, the most relevant excited state for three 
particles is an even doublet. If the approximation is valid, 
the effective matrix element should be smaller than the 
difference between the energy of these states 



43)1 



» M-gl, 



(27) 



where 



We have presented an effective site approach to study 
the conductance through an arrangement of N QD's, di- 
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m: 



off 



j,L'.n=2,4 



(3): 



E, 



(") 



E. 



(3) 



1^, 



(3)\ 



(")\ /,/,(") I 



',(3)\ 



E. 



in) 



E, 



(3) 



(28) 



r 



is the effective matrix element between the ground state 

\tpg ) and an excited state |V'e )• To estimate this matrix 
element in the case of the linear trimer, for simplicity, we 
have restricted the sum above to the lowest states with 
2 and 4 particles in the EHSC. The simplified expression 
is 



/ 1 



E, 



(2) 



E, 



(3) 



E. 



(2) 



E, 



(3) 



with aa,b being the generahzation of Eq. (15) for any 
pair of states \a) and \b). For the parameters used in 

Ei^'^l is 0.16 for 



roff 



\E 



(3) 



Fig. 4, the largest value of Mg 
A/nU = 4.0 X 10-2. 

In the case of the isosceles triangle the above criterion 



is taken at the value of Vg which renders Eg 



(2) 



E. 



(4) 



For the parameters in Fig. 6, the largest value of 



E^^^l is 3.3 X 10-2 for A/ttU = 3.2 x 10" 



M^glJ\Ei'^- 

and for Fig. 8, the corresponding value is 0.29 for 
A/ttU = 1.6 X 10-^. Thus, we conclude that this set of 
parameters is quite consistent with the approximation. 

Quantitative details of our results are also affected by 
two other approximations related with our slave-boson 
treatment: i) the neglect of the triplet states in the sub- 
spaces of A'^ ± 1 particles in the interacting region, and 
ii) the fact that the saddle point chooses the most conve- 
nient solution between those with definite parity in these 
subspaces: either ej+ = dj+ = or ej^ ~ dj^ = for 
all j. However, in the Kondo regime of most interest, 
the model reduces to a one-impurity Kondo model, and 
for the linear trimer, both approximations slightly mod- 
ify the Kondo temperature and in opposite directions. '^^ 
Therefore our conclusions are not affected. In particular, 
the system remains a Fermi liquid with ideal conductance 
in the EHSC. For other problems and larger V, the effects where 



of the triplets can be important near the singlet-triplet 
degeneracy."*^ 

As argued before, if the system is a Fermi liquid, 
its properties at low-energy are expected to be the same 
as those of a non-interacting system, and therefore, one 
expects non-zero conductance in the EHSC. Then, the 
possibility of vanishing conductance in this case, seems 
to be related with the breakdown of the Fermi liquid. 
For a trimer [N = 3) with C^y symmetry, a non-Fermi 
liquid ground state results as a consequence of the de- 
generacy between odd and even doublets in the ground 
state of the interacting region. •^^■■^^ In this case, although 
Eq. (22) suggests a vanishing conductance, our approach 
breaks down because Eq. (27) is not satisfied. However, 
the geometry of the system is incompatible with "trian- 
gular" Csv symmetry, and even a very small deviation of 
this symmetry restores Fermi liquid behavior. 37 there- 
fore the most likely explanation of the result of Biisser 
et al.2* of vanishing conductance in the EHSC for odd 
iV > 1 seems to be a failure of the embedding proce- 
dure for > 2, in spite of its success when only one 
QD is involved. ^^-2^. Recently, the case = 3 has been 
studied by NRG^^ and by Monte Carlo and variational 
techniques. These results agree with ours for small V. 



Appendix A: THE SADDLE-POINT 
APPROXIMATION 



This is a mean field approximation in which the prob- 
lem reduces to minimize the effective free energy F as a 
function of the numbers ej^^,dj^^ and Su- From Section 
lie, F is given by 



-^fermioni 



(Al) 



J 



E 



boson 



E 



+ A') e^] + E [(^^'^ + ^' - ^l] + E {i^S 2i?(^) + A' - 2A.) d^ - A', (A2) 



and 

^fcrmion 



J2 In (l + e-'='(-^))p.(c^) du;, 



(A3) 



where Pa{^) is the total fcrmion density of states given 
by 

Pa{uj) = --ImTrG.H, (A4) 

where the matrix Gcr(w) contains the retarded Green 
functions of the fermionic effective problem. 

It is convenient to separate the fermionic free energy 
into a part i^fermion corresponding to the system without 
the effective site or "impurity" (and therefore indepen- 
dent of the bosonic variables) and the rest, which repre- 



sents the effect of the "impurity" on Fformion 

-^fcrmion ^fcrmion ^ ^-^fcrmion- (-^^) 

The same separation can be made for the density of 
states. Using the relation^'^ 

1 d 

Ap^ = -Im— In G //^ {uj) ( A6) 

the change in free energy after adding the effective site 
becomes 



AFf 



fermion 



(A7) 



where Gffa- is the Green's function of the fcrmions / at the effective site. Since the effective fermionic problem is 
non-interacting, Gffa can be calculated easily using equations of motion. The result is 



Gff'yi^) = lim 



(A8) 



where g2,i/,cr(w) is the unperturbed Green's function 
(without the effective site) at site 2 with parity v and 
spin a. 

Integrating by parts in Eq. (A7) and taking p, ~ 0, 
one obtains 

1 f° 

^ ^ J — oo 



Minimizing the effective free energy, one obtains the 
following equations 



^ 

= 

^ 

= 



dF 

(7 

dP 



Im 



(AlO) 



9do ^ ' 



Vu,aaj,iy duo g2,y,a{'^)Gffa{uj) 



(All) 



dF 

.s.(A'-A.)-^Im|^ 
dF 

E4^ + 2^' + E4.-i> 



K.cr (Q!j,^ej> -I- fij.udj^y) / duj g2.i,M{^)G f f„{uj) 



(A12) 



(A13) 



0,v 



3,v 
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u = 



= r <koGfM-)-sl^J2dl- 



(A14) 



Since one does not expect breaking of SU(2) symmetry 
we take S| = = s and A| = A| = A. For small V we 
can take a constant unperturbed density of states with its 
value at the Fermi energy po — l/(27rt).^^ This simplifies 
the calculations and docs not affect our conclusions. 



g2A^) = -™poO{D - \uj\), 



(A15) 



where 2Dpo = 1 (£) = 7r<). 

For TV = 3, there are 15 independent variables. How- 
ever, with a little algebra the number of independent 
variables can be reduced to four. Using Eq. (AlO), we 
can relate all boson numbers ej^n in terms of one of them. 
We choose eo,,y corresponding to the ground state of iV— 1 
particles: 



(A16) 



where we have defined 



(2) ^ 
'3,1^ 



p(2) 



+ A' 



(A17) 



Similarly, using Eqs. (AlO) y (All) all dj^i, can also 
be related to eo.jy giving 



ai 



En + A' 



2Ei^) - 2A 



where 



R 



(4) 



a 7 



(A18) 



(2) 
0,1/ 



A' 



E. 



(4) 



A' - 2Ei^) - 2A^ 
Using these relations in Eq. (AlO), we can write 



(A19) 



- eo,. <( ( A' + - 2y2s2^^^^^^ 



where (5 = Trpo + Eq. (A14) takes the form 

Finally (A12) becomes 



In 



{E(^)+Xf + 5^ 



tan 

TT 



£;(3) + A' 



= 1 



^Vo(E. [i:,e,,.o.fl (B!fi+B^ 
A' -A 



In 



(i;(3)+A)'+^2 ' 
{E(^) +X + Wf +5^ ^ 



(A20) 



(A21) 



(A22) 



In practice, Eq. (A13) has been used to express s in 
terms of the other bosonic variables and Lagrange mul- 
tipliers. This leads to a system of four coupled non- 
linear equations with the unknowns {eo,+ , eo,-, A, A'}, 
which was solved numerically. In addition, one of the 
eo,i/ should vanish to satisfy both Eqs. (A20). For the 
trimer eo,- = 0. To facilitate the numerical solution, 
we started solving the equations in the EHSC, for which 



Cj^i, = dj^i,, and then A = E^^\ as can be easily seen 
from Eqs. (AlO) y (All). After the solution in this case 
was obtained, the resulting values were used as a starting 
guess for a slightly modified Vg and the process was re- 
peated until the whole range of gate voltages of interest 
was covered. 
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Appendix B: ANALYTICAL RESULTS FOR ONE 
DOT 



As in Rcf. 47; A can be expressed in terms of s and d 
using 



If the interacting system is composed of only one dot, 
in the absence of a magnetic field B, the change in the 
effective free energy in the saddle-point approximation 
depends on two independent variables and the problem 
is simplified considerably. 

For B = 0, S| = S| = s. Eliminating e from the 
constraint + 2s^ + (P = 1 [see Eq. (6)], the change in 
the effective free energy can be written as 



AF 



ID 



ID 



ID 

boson 



boson + ^-^< 
.2 



ID 



fcrmion ' 

2A + [U -2\) (f, 



AF 



ID 

fcrmion 



--Im 

TT 



InGj^A^) dio , (Bl) 



with 



G}'},iLu)^[uj-{Ea + X) + iS]-\ (B2) 
where the half width of the resonance is 

S = As\d+ef. (B3) 



1 r° 1 

s^+d^ = {nf„) = Im / Gif^ = - tan-i 



givmg 



\{s,d) 



tan[7r(s2 + d^)] 



Ea + 6J' 
(B4) 



(B5) 



The two-terminal Landauer formula at T = gives 



G = ^S'\G}'},\l^o = ^ sin' ^W.)), (B6) 



where in the second equality, Eqs. (B2) and (B5) were 
used. Thus, Eq. (26) is verified. 

Replacing Eqs. (B2) and (B5) in (Bl) the change in 
free energy becomes 



AF^^ = Ud^ + 2Ed{s^ + d^) + 



In 



(7^ + 7)2 + ,52 
I 



2(1? + 7) arctan 



D 



(B7) 



where 7 = Scot[T:{s^ + d^)]. In the EHSC, = = 
1/2 — s'^, and the energy depends on only one variable s, 
which can be obtained by minimization. In the Kondo 
limit of large U, d^ 0, the resulting transcendental 
equation can be simplified giving 



which has the correct exponential dependence. If in- 
stead, the mean-field approach of Ref. 40 is used (includ- 
ing some suitable roots to reproduce the non-interacting 
case), an additional factor 4 appears in the denominator 
of the exponent in Eq. (B8). 



D f-TiU 
exp 



4A 



and using Eq. (B3) 



4A 
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